Early time kinetics of systems with spatial symmetry breaking 
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In this paper we present a study of the early stages of unstable state evolution of systems with 
spatial symmetry changes. In contrast to the early time linear theory of unstable evolution described 
by Cahn, Hilliard, and Cook, we develop a generalized theory that predicts two distinct stages of the 
early evolution for symmetry breaking phase transitions. In the first stage the dynamics is dominated 
by symmetry preserving evolution. In the second stage, which shares some characteristics with the 
Cahn-Hilliard-Cook theory, noise driven fluctuations break the symmetry of the initial phase on a 
time scale which is large compared to the first stage for systems with long interaction ranges. To 
test the theory we present the results of numerical simulations of the initial evolution of a long- 
range antiferromagnetic Ising model quenched into an unstable region. We investigate two types of 
symmetry breaking transitions in this system: disorder-to-order and order-to-order transitions. For 
the order-to-order case, the Fourier modes evolve as a linear combination of exponentially growing 
or decaying terms with different time scales. 

PACS numbers; 



I. INTRODUCTION 

Many systems have effective long-range interactions 
and exhibit a variety of crystalline phases. Some ex- 
amples are metals [ij, block copolymers 0, and sys- 
tems with modulated phases d, |3, S]- The kinetics 
of the transitions from a disordered, spatially uniform 
state to a spatially periodic state (disorder-to-order) 
have been explored in several contexts (see, for exam- 
ple, Refs. d, 0, H, Q-) Some aspects of the kinetics from 
one spatially periodic state to another (order-to-order) 
have also been investigated [l^, [ll|, [l^, . 

For systems with long-range interactions the phase or- 
dering process is separated into an early time regime de- 
scribed by a linear theory, an intermediate regime where 
nonlinear effects first become important, and a late stage 
domain coarsening regime. Phase ordering, a process 
whereby a disordered state evolves into a state with 
macroscopic regions of like components, has been well 
studied for systems described by a Landau- Ginzburg free 
energy, such as the ferromagnetic Ising model ■ If the 
coarse-grained order parameter is locally conserved, the 
phase ordering process is called spinodal decomposition. 
If it is not conserved, the process is called continuous or- 
dering. The latter case is referred to as Model A and the 
former is referred to as Model B [l^. 

"Phase ordering" as used here is distinct from a 
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disorder-to-order transition, which we take to mean a 
transition from a disordered to a spatially periodic state. 
A disorder-to-order transition breaks spatial symmetry; 
that is, the initial state has a spatial symmetry that is 
absent in the final state. In contrast, phase ordering for 
the ferromagnetic Ising model with a non-conserved or- 
der parameter has uniform spatial symmetry in both the 
initial and final phases. Thus, in this case phase ordering 
docs not break spatial symmetry [l6l | . 

The early time linear regime of phase ordering was 
first studied by Cahn and Hilliard [O, d, [H and later 
by Cook [iOl. The early time theory was extended to 
the intermediate time regime by Langer, Bar-on, and 
Miller 21] for Model B and by Billotet and Binder ^ 
for Model A. The Lifshitz-Slyozov theory [2^ describes 
the late time regime for Model B which is dominated by 
the surface tension of domain walls. The theory of late 
stage boundary motion for model A was developed by 
Allen and Cahn [2^. 

The Cahn, Hilliard, and Cook (CHC) linear theory 
was originally developed to describe spinodal decompo- 
sition for systems such as binary alloys, where the order 
parameter is locally conserved. The CHC theory pre- 
dicts exponential growth of certain Fourier modes im- 
mediately after a temperature quench into the unstable 
regime [l3, [iM HO] • For systems with long-range in- 
teractions, this exponential growth lasts for a time pro- 
portional to the logarithm of the interaction range |;25|] . 

CHC theory also describes continuous ordering for sys- 
tems such as the ferromagnetic Ising model, with a non- 
conserved order parameter. In this case the temperature 
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is quenched to below the critical temperature at zero ex- 
ternal field (a critical quench). However, if the external 
field is non-zero during the temperature quench (an off- 
critical quench), CHC theory is not applicable. We will 
show that CHC theory also fails to describe a wide class 
of transitions that break spatial symmetry. 

In this paper we extend the CHC theory to describe 
transitions that break spatial symmetry and whose dy- 
namics do not conserve the order parameter (26| . For 
these systems we find two distinct stages of the early 
time evolution in contrast to CHC. Stage one, which is 
absent in CHC theory, is dominated by symmetry pre- 
serving dynamics. The time scale for which this stage 
applies is independent of the interaction range. Stage 
two is analogous to the early time evolution described by 
CHC because the dynamical equations are characterized 
by exponential solutions, and the time scale for which 
stage two is valid is longer for longer interaction ranges. 
For disorder-to-order transitions certain Fourier modes 
grow exponentially. In contrast, for order-to-order tran- 
sitions the Fourier modes evolve as a linear combination 
of exponentially growing or decaying terms with different 
time scales. 

Our particular system of interest is the long-range an- 
tiferromagnetic Ising model. This system has spatial 
symmetry breaking disorder-to-order and order-to-order 
transitions. We use Monte Carlo techniques to simulate 
the system and compare our results to our numerical so- 
lutions of a Langevin equation and to our analytical pre- 
dictions. 

The structure of the paper is as follows. Sectionllllsum- 
marizes the theoretical background on unstable state evo- 
lution for models described by a Landau-Ginzburg free 
energy, focusing on the linear CHC theory. We extend 
this theory to include symmetry breaking transitions in 
Sec. mil The role that symmetry breaking plays in the 
interplay between the two stages is discussed. The long- 
range antiferromagnetic Ising model and its phase dia- 
gram are introduced in Sec. IIVI We develop the gen- 
eralized theory for this model by using the Langevin 
equation based on a coarse-grained free energy as well 
as the Langevin equation describing Glauber dynamics. 
Our numerical and theoretical results are presented for 
the disorder-to-order case in Sec. |V] and for the order- 
to-order case in Sec. I VII Section IVIII summarizes and 
discusses our results. 



II. REVIEW OF CHC THEORY FOR 
CONTINUOUS ORDERING 



In this section we review the Cahn-Hilliard-Cook 
(CHC) theory as it is applied to continuous ordering for 
which the order parameter is not conserved. We begin 
with the Langevin equation for the evolution of the order 



parameter (f)(x,t): 
d(j){x, t) 



-M 



(1) 



where F is the coarse grained free energy, and M is the 
mobility (which for the CHC theory is assumed to be a 
constant in space and time) and 



B = 2MT 



(2) 



from the fluctuation dissipation relationship where we 
have set ks — Equation ([1]) is also known as the time- 
dependent Ginzburg-Landau equation The first 
term on the right-hand side of Eq. ([1]), which we call the 
drift term, lowers the free energy as time evolves. The 
second term describes the noise and represents the small 
scale fast dynamics as random forces. We will take the 
noise to be white Gaussian and to obey the conditions: 



{t]{x, t)r]{x', t')) = 6{x - x')6{t - t'). 
By using the the Landau-Ginzburg free energy, 

F[^{x)]^jdx[^[Vi{x)Y+'-^H^) 
1 



(3) 
(4) 



+ -^\x)-h4>{x)), 



(5) 



Eq. H]) becomes, 
d(j){x, t) 



-M{^R'^W'^^{x, t) + e4>{x, t) 
+ i>^{x,t)-h} + y/BT]{x,t), (6) 



where e — {T ~ Tc)/Tc is the reduced temperature, Tc 
is the critical temperature, h is the external field, and R 
is the interaction range. The extrema of the free energy 
can be found by solving the Euler-Lagrange equation: 
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i?^V>(a;, t) + £(j){x, t) + <j)'^{x, t) 



(7) 



For a critical quench the initial state is equilibrated 
at initial temperature > in zero external field 
/li = 0. The temperature is instantaneously quenched 
to a final temperature Tf < Tc in the unstable regime 
where £/ = {Tf — T^jT^ < 0, keeping the external field 
hf — 0. Before the quench, the system is in a disordered 
state corresponding to the average order parameter, i^q = 
{4>{x)) = for spatially homogeneous systems. We will 
consider small fluctuations, 'ip{x,t), of the order parame- 
ter about this value, (f>{x,t) — (t>Q^ip{x,t) = t/j{x,t), and 
keep terms flrst order in "0(0;, t) so that Eq. ([6]) becomes: 

d^{x, t) 



dt 



-M [-R^V'^'i;{x,t) - \ef\i^{x,ty^ 
+ ^/B'n{x,t). (8) 
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We have used the fact that (?!)o = is a solution to Eq. ([7]) . 
In Fourier space Eq. ([5]) becomes 



dtp{k,t) 



dt 



with 



i^chc(fc)--M(i?^fc^-|e/|) 



(9) 



(10) 



We see from Eq. ^ that the individual Fourier modes 
evolve independently of one another. The solution to 
Eq. dS]) is 

iP{k,t) = V;(fe,0)exp(i:'chc(A:)t) 

+ Vb[ Tj{k,t)G{t~t')dt\ (11) 
Jq 

where G{t - t') = exp((i - t')Dchc{k)). We use the 
definition of the equal time structure factor S{k,t) = 
{(t>{k,t)4){—k,t)) /V, Eq. and the properties of the 
noise (Eqs. ([3]) and ([4])) to obtain: 

S,M) = e-.«'(s,M^O, + -^) 



M 



(12) 



where (if ^TJ^. 

Systems with i? ^ 1 are said to be near- mean-field and 
the limit i? ^ oo is the mean- field limit |28i. For Tf < 
Tc simulations have shown that Eq. is applicable 
for systems with i? ^ 1 for some initial time after the 
quench [H, H HI, lH . Binder pl has argued that the 
linear theory should break down when the fluctuations 
are large in comparison to the average order parameter. 
He predicted that CHC theory breaks down at a time r 
that scales as r ~ In/i (see Eq. ([T^ V 

To understand Binder's result, we consider how V'(fc) 
depends explicitly on R, the relevant length scale of 
the problem. We scale all lengths with R by defining 
r = x/R, and q = kR. We see from Eq. (|10l) that £'ciic(<z) 
is independent of R. It follows that G{t — t') is also in- 
dependent of _R, and from Eq. PT|) ip{k,t), ?/j(fc, 0), and 
?7(fe, t) must depend explicitly on R in the same way. 
Equations ^ and (U) become 



(77(r,t)> 
{r]{r,t)ri{r',t')) 





6{r 



^')5{t--t'). 



(13) 
(14) 



From Eq. (I14p and 5{x) — R '^6{r) the noise must satisfy 
T]{x,t) — R^'^/^r]{r,t), or 



l{k,t) = R-'/^7j{q,t), 



(15) 



where r]{q,t) ^ J dxc ''"'^ 77(r, t). This motivates us to 
introduce ipiq^t) through 



and il){r,t) though ip{q,t) = / dxe~^'''^ ip{r,t) and is 
independent of R. We consider the dynamics oiip{r). If 
we substitute Eqs. ([TS]) and (fTBjl into Eq. ^ and use a 
Laplacian with respect to scaled coordinates = i?^V^ 
all R dependence in Eq. ([8]) cancels. Thus the solution to 
Eq. ([5]), i/j{r,t), is also independent of R. Equations (fT5|) 
and mean that larger R results in effectively smaller 
noise and noise induced fluctuations. 

We may drop the i/j^{x) term in Eq. ([8]) only if iIj{x) ^ 
ip^^x). Thus we estimate the breakdown time r from the 
condition ip'^ (x) ~ 1. For a given value of q this condition 
gives 



which yields 



exp (2J(g)r) 



dlni? 
2D{q)- 



(17) 



(18) 



where we have used Eq. (jT6| along with the time depen- 
dence of Eq. pT|) . Note that the result for the breakdown 
of the linear theory is a consistency argument, and it is 
possible for the linear theory to fail earlier than the time 
given in Eq. (fTS]) . 

The result that the duration of applicability of the lin- 
ear theory is longer for larger R is illuminated by a sim- 
ple mechanical analogy (see Fig. [T|). The CHC theory 
describes the evolution of certain fluctuations of the sys- 
tem about an unstable stationary point of the free en- 
ergy. The system is analogous to a ball sitting at the 
top of a fairly flat hill subject to small random pertur- 
bations. These perturbations will eventually cause the 
ball to roll off the hill; the smaller the perturbations, the 
longer it takes for the ball to roll off. As R increases the 
noise effectively decreases and the system remains near 
the unstable point of the free energy longer. When the 
free energy can no longer be approximated quadratically, 
the CHC theory is no longer valid. 



III. GENERALIZATION OF CHC THEORY 

Equation ([T]) is an example of a dissipative Langevin 
equation. Because we will also consider Glauber dynam- 
ics, we consider a Langevin equation of the more general 
form, 



90(r,^) 
dt 



e[0, r] + VS[0,r]i?-'^/2,7(r, t), (19) 



V'(fe,i) = i?-'^/V(g,t). 



(16) 



where (j){r,t) is the order parameter with scaled coordi- 
nates, the drift term 0[(/),r] is nonlinear in (j){r,t) and 
the noise satisfies Eqs. and p^ . Equation p^ is 
written in scaled coordinates so that the i?-dependence 
of the noise term is explicit (8[0, r] and r] are in- 
dependent of R.) 

In Sec. [IT] we derived the CHC theory as an expan- 
sion about a uniform order parameter background. This 
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background configuration was the solution to the Euler- 
Lagrange equation for the parameters both before and 
after the quench and is therefore time independent after 
a critical quench. More generally, we will expand about a 
background configuration 0b (r, t) which is not necessarily 
uniform and may be time dependent after the quench. 

At the time of the quench, the system is in a con- 
figuration consistent with the equilibrium state before 
the quench. We can write this state as (j){r,t = 0) = 
(j)b{r, t = 0) + R-'^/'^ip{r, t = 0). The background config- 
uration (l)b{r,t = 0) is a solution to the Euler-Lagrange 
equation for Ti and hi, is independent of the noise, and 
may be nonuniform. The term R~'''^^ip{r,t = 0) repre- 
sents small fluctuations about the background due to the 
noise. In general, (t>b{r,t = 0) 7^ </'o('"); that is, the back- 
ground configuration is not a stationary point of the free 
energy and may be time dependent after the quench. We 
are interested in how this background configuration and 
these fluctuations evolve with time after the quench, that 
is the time evolution of 



(r,t) =0b(r,t) + i?-''/V(r,<). 



(20) 



If we expand Eq. to first order in the small fiuctua- 
tions, R~'^^'^tlj{r,t), about the background configuration 
(j)b{r,t), we obtain 



dcj)b{r,t) 



dt 



dr 



dt 

,S&[4,,r 



54>{r\t) 



i?-'^/2i/;(r',t) 



/B[0~r] 



4>=4'b 
R-'^/^r]{r,t). 



(21) 



It is natural to separate Eq. ipT|) terms of order from 
terms of order R~'^/'^ giving the dynamical equations for 
the background and the fluctuations. The equation for 
the evolution of the background configuration, which is 
zeroth order in i?^*^/^, is 



dt 



= Q[(j)]\4,=4,h- (evolution of background) (22) 



Equation (|22|) is noiseless and independent of R. In con- 
trast, the equation for the evolution of the fluctuations 
is noise driven and is first order in R^'^^^: 



d^{r,t) 
dt 



dr 



i>{r\t) 



ir',t) 



(evolution of fluctuations) (23) 

The background configuration 0f,(r,t) and the fluctua- 
tions ip{r,t) evolve simultaneously; Eqs. ([H]) and ([^5]) 
together describe the early time kinetics of unstable sys- 
tems. Equation l|23p differs from the CHC theory de- 
scribed in Sec, mi because the coefficient of ip{r,t) on the 
right-hand side of Eq. depends on (t)b{r,t), which 
evolves according to Eq. (1^^ . Because of the explicit 



time dependence of (j)b{r,t), we do not expect to find 
exponential solutions as found in the CHC theory. 

Equation ((22)) dominates the evolution of (j>{r,t) until 
1 6 [0] 1 0=0 J ~ i?"''/^ (we wifl discuss specific examples of 
this evolution in the following). We refer to the kinetics 
while the background configuration is evolving as stage 
one. If 



|e[<^]| 



0=06 



< R 



-d/2 



(24) 



we refer to the kinetics as stage two. 

If 8[0] cx SF/Scf), then Eq. becomes a dissipative 
Langevin equation. In the limit i? — s- 00 the condition for 
stage two evolution (Eq. (j24p ) implies that SF/Scf) ~ 0, 



Mr,^) = '/>o(r). 



(25) 



where (po{r) is an unstable solution to the Euler-Lagrange 
equation for the final parameters Tf and hf corre- 
sponding to the temperature and external field after the 
quench, respectively. In stage two the dynamics shares 
several features with the dynamics of CHC. Namely, the 
background configuration is nearly a stationary point of 
the free energy, and the linear evolution admits (nearly) 
exponential solutions, which are applicable for a time 
T Ini?, as in Eq. (fTO|) . 

If 4)b{r,t = 0) ^ 0o(^), then stage one kinetics begins 
immediately after the quench, and either there is no stage 
two kinetics at any time after the quench, or stage two 
kinetics follows stage one. If (j)b{r,t = 0) = (f>Q{r), then 
stage one kinetics is absent and stage two kinetics begins 
immediately after the quench. Table U summarizes the 
various possibilities. 

A critical quench in the ferromagnetic Ising model (see 
Table H] and Fig. 1(a)) exhibits stage two but not stage 
one kinetics. In this case Ti — 00, hi = 0, Tf < Tc, 
hf —0, and 0f,(t = 0) = 0o- 

In contrast, an off-critical quench in the ferromagnetic 
Ising model (see Table |I] and Fig. 1 (b) I is an example 
that exhibits stage one but not stage two kinetics. In the 
off-critical quench the initial conditions are the same as 
for the critical quench [Ti = 00 and hi = Q), but the sys- 
tem is quenched below T^ and hf ^ 0. Because hf ^ 0, 
the initial background configuration (f>b(t = 0) = is not 
a solution to the Euler-Lagrange equation for the final 
parameters Tf and hf. The background configuration 
evolves via Eq. (|22p directly toward the stable minimum 
(f)m of the free energy landscape and away from the un- 
stable solution (po (see Fig. 1(b) ). Hence, the system will 



never reach the unstable stationary point, and there will 
be no stage two kinetics. 

In contrast to the ferromagnetic Ising model [s^, we 
find stage one followed by stage two kinetics in systems 
that break spatial symmetry. In this case the noiseless 
dynamics in stage one and the initial configuration con- 
spire to drive the configuration toward the unstable so- 
lution of the Euler-Lagrange equation 0o- In stage one 
both (t)b{r) and ip{r,t) evolve simultaneously until the 



5 



Case 


Spatial symmetry breaking 


Stage one 


Stage two 


Example 


Mt = 0) ~ 00 


No 


No 


Yes 


Critical quench for ferromagnet 


Yes 


No 


Yes 


Critical quench for antiferromagnet 


Mi = 0) / 00 


No 


Yes 


No 


Off-critical quench for ferromagnet 


Yes 


Yes 


Yes 


Off-critical and external field 
quench for antiferromagnet 



TABLE I: Summary of early kinetics for various types of transitions for the long-range ferromagnetic Ising model and the 
long-range antiferromagnetic Ising model. In the table 0(,(t — 0) is the background configuration just after the quench and 0o 
is the unstable solution to the Euler-Lagrange equation for T/ and hi. These quenches are defined in the text. 




(b) 



FIG. 1: A representation of the free energy landscape for a 
ferromagnetic Ising model as a function of the average mag- 
netization (0) for (a) a critical quench and (b) an off-critical 
quench for hf > 0. The dot represents the background con- 
figuration of the system immediately after the quench. For 
a critical quench the system is initially at an unstable sta- 
tionary point of the free energy, and the CHC theory governs 
its early evolution. For long-range systems the noise is effec- 
tively small, and the system fluctuates about the top of the 
hill for an extended time. Eventually, the system will be in 
one of the two global minima (j>m if the system size and tem- 
perature are finite. For an off-critical quench the background 
configuration sits to the right of the free energy stationary 
point. The noiseless dynamics of Eq. (|22p drives the system 
directly toward the stable minimum 0m and away from the 
stationary point. Thus, in this case we do not expect stage 
two exponential growth at any time after the quench. 



system reaches the unstable solution 0o for Tf and hf. 
Stage one kinetics is analogous to a ball rolling down a 
hill toward a saddle point (see Fig. ^ . When the system 
reaches the stationary point, stage two begins. Stage two 
is analogous to the ball fluctuating about the top of the 




long 
ranges 



FIG. 2: A representation of the free energy landscape for 
a system with spatial symmetry breaking. For systems with 
short range interactions fluctuations cause the ball (represent- 
ing the background configuration of the system) to roll off the 
hill, the ball does not approach the saddle point and stage two 
behavior is not observed. In contrast, for systems with long- 
range interactions the evolution of fluctuations governed by 
Eq. (|23|l is suppressed by a factor of i?"'*''^, and the noiseless 
dynamics of Eq. (|22|l dominates. Due to this noiseless sym- 
metry preserving dynamics the ball rolls down the hill toward 
the saddle point along the dashed line with little deviation 
due to the fluctuations. The dashed line represents all con- 
figurations with the same symmetry as the initial state. The 
directions perpendicular to the dashed line represent symme- 
try breaking degrees of freedom. Once near the saddle point, 
the ball essentially stops rolling (condition (|24p is satisfied) 
and stage two kinetics takes over as unstable fluctuations grow 
for a time of order r ~ Ini?. Eventually, after stage two, the 
ball will roll off the hill along the dotted line. 



hill and ends when the ball finally rolls off the top of the 
hill towards a stable state. 

Why does the ball in Fig. [2] not roll off the hill be- 
fore it reaches the saddle point? Symmetry considera- 
tions provide the answer. For dissipative dynamics where 
Q oc -~6F/5(j), the noiseless dynamics of Eq. (|22|) pre- 
serves certain spatial symmetries, depending on the form 
of the free energy. The dynamics of Eq. ([22|) evolves 
06(r,i) to the lowest free energy within these symmetry 
constraints. This configuration is an unstable stationary 
point of the free energy. It is a local minimum with re- 
spect to the symmetry constrained degrees of freedom, 
but a local maximum with respect to the other symme- 
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try breaking degrees of freedom. Once at this saddle 
point, the nonhnear background dynamics of stage one 
terminates, and stage two kinetics takes over. See also 
Ref. [H. 

The spatial symmetries preserved by the noiseless dy- 
namics of Eq. (I22|) arc not preserved by the noise driven 
dynamics of Eq. that is, the noise causes ran- 

dom fluctuations which break the spatial symmetries of 
06 (r, t). Initially, these spatial symmetries are broken by 
'4'{r,t = 0), the fluctuations about the background at 
the time of the quench due to ri{r,t < 0), the noise in 
the system before the quench. These initial fluctuations, 
which are of order R~'^/^ , begin to grow (or decay) imme- 
diately after the quench according to Eq. (^5)) . However, 
during stage one the fluctuations ^-!{r, t) remain of order 
R-^/^. The reason the ball in Fig. [2] does not roll off the 
hill is because the fluctuations which would drive the ball 
off the hill are small for sufficiently large R. The time it 
takes for the ball to roll down to the saddle point (the 
duration of stage one), is independent of R. In stage two 
kinetics, when the nonlinear dynamics of Eq. (|^^ has 
ceased, the symmetry breaking fluctuations continue to 
grow for an extended period of time, eventually becoming 
of the order 1. 



Clumps 




0.05 0.1 0.15 0.2 
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FIG. 3: Approximate phase diagram for the long-range anti- 
ferromagnetic Ising model with a square shaped interaction. 
The system has two low temperature ordered phases, the 
clump phase and the stripe phase, and a disordered phase. 
The darkened square in the snapshot of the disordered Ising 
model shows the interaction range for a spin at the center of 
the square. 



IV. THE LONG-RANGE 
ANTIFERROMAGNETIC ISING MODEL 

A. The model and phase diagram 

To investigate a system with spatial symmetry break- 
ing phase transitions, we simulate the long-range anti- 
ferromagnetic Ising model using Glauber Monte Carlo 
dynamics [36| . At each update a spin is selected at ran- 
dom and flipped with the transition probability W = 
(1 -|-e'^^^)"\ where l^E is the change in energy incurred 
from flipping the spin. Time is measured in units of 
Monte Carlo step per spin. For N spins in a system, 
one Monte Carlo step (MCS) is equal to N spin flip at- 
tempts. With this definition, it is natural to consider 
fractional Monte Carlo steps in which fewer than N spin 
flips have been attempted. 

We study the two-dimensional long-range antiferro- 
magnetic Ising model with a square shaped interaction. 
For this interaction a spin at a; = {x,y) interacts with 
a second spin at x' — {x',y') only if \x — x'\ < R and 
\y-y'\<R{see Fig.©. 

There are two types of periodic low temperature con- 
flgurations, which we call the stripe phase and the clump 
phase, that are energetically favorable with respect to the 
disordered phase. The system can lower its overall inter- 
action energy in the clump or stripe phases by grouping 
like spins together. Although a spin pays an energy cost 
due to its proximity to the other spins in its own clump 
or stripe, it avoids interacting with most of the spins in 
neighboring clumps or stripes. This system is related to 
a previously studied model of particles [IJ- 



An approximate phase diagram is shown in Fig. [31 
The pictures are snapshots of the conflgurations from 
our Monte Carlo simulations. From Fig. [3] we see that 
the disorder-to-clumps transition as well as the disorder- 
to-stripes transition are symmetry breaking disorder-to- 
order transitions. In both cases the uniform spatial sym- 
metry is lost. The two order-to-order transitions are also 
symmetry breaking. In the transition from stripes to 
clumps the continuous translational symmetry along the 
direction of the stripes is lost. In the clumps to stripes 
transition the 90° rotational symmetry is lost. 

The orientation of the square interaction restricts the 
possible alignments of the stripes and clumps: stripes 
must be either horizontal or vertical and clumps must 
be aligned in the vertical and horizontal directions. A 
circular-shaped potential which we do not consider here 
gives a honeycomb arrangement of clumps and stripes 
with no angular restrictions. 

The Hamiltonian for the long-range antiferromagnetic 
Ising model is 

= -— ^ CTiCTj - ft-^CTi, (26) 

<ij> i 

where the brackets indicate that the summation is over all 
spins within the interaction range R, and ai = ±1 refers 
to the value of a spin at site i. The interaction J = ±1, 
where the positive sign corresponds to the ferromagnet 
and the negative sign corresponds to the antiferromagnet; 

z = (2i? + l)'' - 1 (27) 
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is the number of spins within the interaction region of 
a single spin. ScaHng the strength of the interaction by 
the number of the spins in an interaction range keeps 
the interaction energy per spin finite as we increase the 
R 



which we are free to choose to correspond to the par- 
ticular system of interest. We allow the mobility to be 
spatially varying. We take 



ei 



-M(r) 



SF 



(31) 



B. Dissipative Langevin dynamics 

We now develop the equations for the relaxational 
Langevin dynamics as well as the background and fluc- 
tuation dynamics for a free energy which corresponds to 
the long-range antiferromagnetic Ising model. We begin 
with the coarse-grained free energy 



F 



[^(r)]=i?''(-^y dr((A*</.)(r))0(r) 

dr{Ts[cl,{r)]+hcP{r))), (28) 



where 



1 



r = 



(A) 



2 Si is the coarse-grained magne- 



tization of the spins inside a box of length A in c? di- 
mensions centered at position r. We use the notation 
(A * (/))(r) = J dr'A{r — r')(t){r') to denote a spatial con- 
volution. As in Eq. ([26|l J = +1 for the ferromagnet and 
J = — 1 for the antiferromagnet. The entropy density 
s[0(r)] due to the coarse-graining procedure is 

1 



s[0(r)] = 21n2--{(l 



(r))ln(l + (/)(r)) 
+ {l-cj,{r))\n{\-(t>(r))}. (29) 



The scaling of the strength of the interaction is implicit 
in the interaction function A(r) which we take to be a 
repulsive pairwise step potential of the form 



{-R<x,y<R) 
otherwise 



(30) 



where z - {2R)'^ for large R (from Eq. This Kac 

form of the potential [38| | yields a Hamiltonian that is 
well defined in the mean field limit; the energy per spin 
remains finite as i? ^ oo because / drA{r) is independent 
of R. We believe that this form of the free energy exactly 
represents the Ising model because a similar expression 
gives the exact free energy for fluid density derived by 
Grewe and Klein [Uli^. 

We can recover the Landau-Ginzburg free energy func- 
tional in Eq. ([5]) by expanding A(r — r') about the q = 
mode and neglecting terms higher than 0{q^), and also 
expanding the entropy term in powers of (j) and truncat- 
ing terms higher than 0((/)^). For the antiferromagnet the 
Fourier modes of interest are not centered about q = 0, 
so we will not truncate the potential. 

For dissipative Langevin dynamics the drift term 
0[(/),r] is proportional to —SF/6(f> times a mobility M, 



where F = F[(f>{r)]/ R"^ is the free energy per interac- 
tion volume for large R and is independent of R since 
F[(t){r)] (X R'^ (Eq. The equation of motion, 

Eq. (fTO|) . becomes: 

^^1^ = M(r){ J(A * 0)(r,t) - arctanh[0(r, t)] 

+ h^ + ^B{^R-'^/^rjir,t), 

(dissipative Langevin dynamics) (32) 

where B{r) = 2Af(r)/3^^, as determined by a fluctua- 
tion dissipation theorem. We numerically solve Eq. (j32|) 
with simulated noise in order to compare the relaxational 
Langevin dynamics with our predictions for stage two ki- 
netics. 

The dynamics in Eq. (|32p evolves the system toward 
a local free energy minimum which is a solution to the 
Euler-Lagrange equation SF/S(p = 



0o(r-) - tanh [/3J(A * ?io)(r) + . 



(33) 



To compare to the Glauber dynamics in Sec. IIVC| we 
now expand about the background configuration to find 
equations for the evolution of the background configura- 
tion and the noise driven fiuctuations. Using the form 
of 8[0, r] given in Eq. PT|) . the nonlinear background 
dynamics in Eq. ((22l) becomes 



. M(r)(j(A*0,)(r-,t) 

- PJ^ arctanh[(/)f,(r, t)] + , (34) 

and the linear theory for the fluctuations represented by 
Eq. (^5]) becomes 



dt 



M(r.)(j(A.^)(r,.)-/?-^^|^ 



yjB{r)r){r,t). 



C. Glauber Langevin dynamics 



(35) 



There is a Langevin equation corresponding to the long 
range Glauber Monte Carlo dynamics which is not equiv- 
alent to the dissipative Langevin dynamics of Ecu (I32p. 
This Glauber Langevin equation, derived in Ref. [41|, is 



dt 



tanh [J3J{K * 0)(r,i) + f3h 



(36) 
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where B[<j)] = 2 - tanh2(/3J(A *(j)){r,t)+ ph) - (f>'^{r, t) 
depends explicitly on (f){r, t) and the noise obeys the con- 
ditions in Eqs. ITS]) and (HI]). See also Ref. [ia ]. 

Glauber dynamics obeys detailed balance and there- 
fore should take the system to equilibrium. We see that 
if we set the noise part of the dynamics equal to zero 
the system converges to a solution of the Euler-Lagrange 
equation 



= tanh(/3J(A * cj)o){r) + (ih) - Mr), (37) 



dt 



in agreement with Eq. (j33p . 

Following a quench to (3f and hf the linearized equa- 
tion for Glauber dynamics (Eq. ([36|l ) from Eq. ([23]) is 

X Jl3f{A*tlj){r,t) -'ip{r,t) 



- yS|^I^?7(r,i). 



(38) 



When (j)}j{r,t) ~ 0o('") we can apply Eq. ([57)) to Eq. 
and obtain an expression for stage two kinetics 



dt 



{l~4>l{r)) f3fJ{A*^){r,t)~i^{r,t) 



(39) 



Equation (j39p for Glauber Langevin stage two dynam- 
ics is identical to Eq. (f35|) for the dissipative Langevin 
dynamics of stage two if the mobility is chosen to be 



M{r)^Pf{\-ct>l{r)). 



(40) 



V. DISORDER-TO-ORDER TRANSITIONS 
A. Theory 

We first treat the transition from the disordered phase 
to either stripes or clumps in the antiferromagnetic Ising 
model. We consider dissipative dynamics (Eq. ([M]) ) 
where we choose the mobility to obey Eq. (|40)) . where 
(/)o(r) = 00 so that M{r) = M . For this choice, the dis- 
sipative dynamics and the Glauber dynamics should be 
equivalent for stage two. 

For comparison we will also discuss the contrasting be- 
havior of the ferromagnetic Ising model. In this case the 
initial background configuration is uniform and 0f, {t = 0) 
may or may not be a solution to the Euler-Lagrange 
equation. Consider the uniform solutions to the Euler- 
Lagrange equation (Eq. (I55|) ) given by 



00 = tanh(/3(J0o + h)) , 



(41) 



where we have used J drh.{r) — 1. For ft, = 0, 0o = 
is a stable minimum for the antiferromagnet and the 
ferromagnet at high temperatures. For the ferromagnet 



at low temperatures and relatively small there is one 
solution unstable to g = fluctuations and two solutions, 
either both stable, or one stable and one metastable. 

For the antiferromagnet at low temperatures and rel- 
atively small we know that the stable solution is not 
spatially uniform. The uniform unstable solution given 
by Eq. is stable for the q = mode and unstable for 
certain > modes. The details of the instability are 
discussed in the following. 

The linearized equation for stage one kinetics is of the 
form of Eq. ([55]) . with (f>b{r,t) = (j)b{t). Equation 
can be written in Fourier space as 



where 



at 



D{q,t) 



D{q,t)^j;{q,t) + VB7^iq,t), (42) 



M {^JA{q) 



(43) 



The solution to Eq. ((32) is not an exponential in general, 
but can be solved numerically by simultaneously solving 
Eq. dHl for Mt)- 

Eventually, we expect that 4>b{t) ~ 0o and that stage 
two kinetics takes over. The dynamics of the fluctuations 
is then given by Eq. ([1^ with (f>i,{r,t) = (f)Q, so that 
D{q,t) for stage two kinetics becomes 



D{q) = f3fil-<j,l)JA{q)-l. 



(44) 



With these simplifications, we solve Eq. (|42)) to obtain 
a generalization of Eq. (fT^ for disorder-to-order transi- 
tions 



S{q,t) = 



5{qf 



V 
1 - 



,2D(q)t 



^0(9) 



l-0j 



D{q) 



D{q) 



(45) 



where So{q) — (\ip{t — to)P) /V is determined by the ini- 
tial conditions for the differential equation (|42)) and can 
be treated as a fitting parameter unless stage two kinet- 
ics begins immediately after the quench where So{q) = 
S{q,t = 0). Equation ([l5|) is only valid for stage two 
kinetics, that is, for t > to. 

The sign of D{q) in Eq. (|44p determines whether the 
structure function grows or decays exponentially. To see 
for which q value D{q) first changes from negative to 
positive as the temperature is lowered to less than Tc, we 
look at the potential A in Fourier space. For the square 
shaped interaction we have 



sm{qx) sin((?y) 



qx 



(46) 



We plot the function D{qx,0) in Fig. g) 

For the ferromagnet (J = -1-1) D(q) for q = first 
crosses from negative to positive as the temperature is 
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FIG. 4: The function D{qx,Q) is plotted for (a) the ferro- 
magnet (J = 1) and the (b) antiferromagnet (J = —1) for 
various temperatures. Using Eqs. (|44p and (I46p the Fourier 
modes corresponding to positive values of D{q) are exponen- 
tially growing. As the temperature is lowered, the ferromag- 
net first becomes unstable at g = 0, but the antiferromagnetic 
becomes unstable at g = g* ~ 4.4934 (see Eq. ^^). 



lowered. As the temperature is lowered for the antiferro- 
magnetic case [J — —l), D[q) first crosses from negative 
to positive at the global minimum of A(q) which occurs 
at q = (±g*,0) or q = (0,±9*) where q* ~ 4.4934. The 
temperature at which this crossover occurs is the critical 
temperature Tc for a given (/jq- The critical temperature 
for the antiferromagnetic system is 



sm{q* 



(47) 



Below this temperature there is at least one unstable 
mode which grows exponentially. 



We find from numerical solutions of Eq. ([32| that the 
disorder-to-order and order-to-disorder transitions occur 
at the same temperature, indicating the existence of a 
continuous phase transition. Thus the line in the phase 
diagram (Fig. [3]) separating the disordered phase from 
the ordered phases is a line of critical points. 



MC data 
DL data 
Theory 




FIG. 5: Structure factor S{q, t) = [S{±q, 0; t) + S{0, ±q; t)]/4 
versus q for the antiferromagnetic Ising model for early time 
t — 0.125 MCS after a critical quench from Ti = oo tp 
Tf = 4rc/9. The Monte Carlo data (MC data) has an in- 
teraction range of _R = 128, while the numerical solutions to 
the dissipative Langevin Eq. p2|l (DL data) has a range of 
lO''. The theory line is a plot of Eq. (|45)l. 



B. Simulation Results 

1. Critical Quench 

We first consider the critical quench for the antiferro- 
magnet for which we expect stage two kinetics immedi- 
ately following the quench (refer to TableU]). For this case 
M = (3f from Eq. (140]) . Figure [5] shows the correspon- 
dence between the linear theory prediction Eq. (P5|) . the 
numerical solutions to the dissipative Langevin equation, 
Eq. ([5^ . averaged over the noise, and the Monte Carlo 
simulations averaged over the noise. In Fig. [5] the struc- 
ture factor is plotted versus q at an early time t = 0.125 
MCS (Langevin time is measured in units of MCS) . Be- 
cause we use a square interaction, we plot 



Siq,t) = [S{±q,0;t) + S{0,±q;t)]/4 



(48) 



instead of taking a circular average. For all simulations 
presented, we hold the ratio of the system size to the in- 
teraction range constant at L/R ~ 2.78 so that the stripe 
state of the final phase is approximately commensurate 
with the system size j43| . 

The coarse-graining employed in the Langevin equa- 
tion allows us to consider system sizes and interaction 
ranges which are much larger than those of the Monte 
Carlo simulations. In Fig.[5]the Monte Carlo simulations 
have R = 128, and the Langevin solutions have R ~ 10^. 
The theory curve corresponds to an infinite interaction 
range (so that the linear theory lasts forever.) Figure [H] 
shows that the linear theory is still a good approximation 
at time t = 0.125 for the Langevin and the Monte Carlo 
dynamics. 

Next we compare the time dependence of S{q, t) given 
by Eq. (|45)) following a critical quench to our Monte 
Carlo simulations. Figure [H] shows S{q, t) for q ~ q* 
plotted against time for several values of R along with 
the theoretical prediction of Eq. (^5)) . As predicted by 
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FIG. 6: Semi-log plot of the structure factor S{q, t) for 
q — 4.516 ~ q* versus time for our Monte Carlo simulations 
of a critical quench in the antiferromagnetic Ising model. The 
temperature is quenched from Ti — oo to Tf — 4Tc/9. The 
theory line is a plot of Eq. (|45p . The duration of the expo- 
nential growth increases as the interaction range increases. 



Binder [25|, we find that stage two kinetics holds for 
longer times as we increase the interaction range for this 
value of q. 



2. Off-critical quenches 

For external field quenches or off-critical temperature 
quenches we find that the short time behavior is indepen- 
dent of R, and corresponds to our prediction of stage one 
kinetics. In Fig. [7| we present the results of off-critical 
quenches {T, = oo, h, ^ 0, Tf = 4Tc/9, hj = 0.5) 
for R = 23, 92, 368 for Monte Carlo simulations av- 
eraged over at least 700 independent initial configura- 
tions for each value of R. We plot S{q, t) versus time for 
q = 4.516 ~ q*. We interpret the first Monte Carlo step 
to be stage one kinetics, which, as we see from Fig.[7]has 
a time scale that is independent of the interaction range. 
At about one Monte Carlo step stage two kinetics take 
over, and the structure factor at g ~ g* grows exponen- 
tially until stage two kinetics breaks down. The break 
down occurs later for longer interaction ranges. These 
observations are consistent with the prediction that the 
time scale of stage one kinetics is determined by the noise- 
less dynamics of Eq. ([221) which is independent of R. Note 
that stage one kinetics were not predicted by the CHC 
theory. 

To show that our generalized theory predicts the same 
behavior for both the disorder to stripes and the dis- 
order to clumps transitions, we compare S±{q,t) (q in 
the direction perpendicular to the final stripe phase) 
and (g, t) {q in the direction parallel to the final 
stripe phase) where S±{q,t) — S{q*,0;t) and S\\{q,t) — 
S{0,q*;t) for vertical stripes and S±{q,t) = S{0,q*;t) 
and S'li = S{q*,0;t) for horizontal stripes. We find that 
S±{q,t) and S\\{q,t) both grow exponentially with the 
same rate during stage two (see Fig 8(a) ). It is not until 



FIG. 7: Semi-log plot of the structure factor S{q, t) for 
q = 4.516 ~ q* versus time for Monte Carlo simulations of 
an off-critical quench in the antiferromagnetic Ising model. 
The simulation parameters are Ti = oo, T/ = ^Tc, hi = 0, 
and hf = 0.5. The first Monte Carlo step is stage one kinet- 
ics, whose time scale is independent of interaction range. The 
duration of stage two kinetics, marked by exponential evolu- 
tion, increases as the interaction range increases. The theory 
line is a plot of Eq. (|45p (with S(i{q) fit to data) and is only 
expected to be applicable for stage two. 



after stage two that S'li (g, t) begins to decrease toward its 
final value. For comparison, S{q, t) is plotted in Fig. |8(b)| 
for the disorder to clumps case. Note that the qualita- 
tive behavior for the disorder to stripe and the disorder 
to clump transitions are the same. 



VI. STRIPES TO CLUMPS TRANSITION 

We now consider the unstable kinetics for spatial sym- 
metry breaking order-to-order transitions. In this section 
we develop a theory that describes the early time kinet- 
ics from an initial stripe state to the stable clump phase 
following an instantaneous quench of the external field. 
The development of the theory is analogous to that of 
the disorder-to-order case, but the periodic structure of 
the initial configuration of the system introduces some 
interesting deviations from the CHC theory. 



A. Theory 

After an external field quench with the temperature 
held constant, the stable solution to the Euler-Lagrange 
equation is a clump configuration corresponding to Tf 
and hf (see Fig.|31). We find from numerical solutions to 
Eq. ((32)) that there is a second, unstable stripe solution 
as well, (f>s{f) = (f'sif-x) where and ry represent scaled 
spatial coordinates (we assume here that the stripes are 
in the vertical direction). 

As in Sec. El we proceed by considering dissipative dy- 
namics (Eq. (|36p ) where we choose the mobility to obey 
Eq. (|40p. where we now choose (f>o{r) = (psirx) so that 
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Time (MCS) 
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FIG. 8: (a.) A semi-log plot of S±{q,t) (where q is in the 
direction perpendicular to the final stripe phase) and {q, t) 
(where q is in the direction parallel to the final stripe phase) 
versus time for an off-critical quench from disorder to stripes 
with Tf = 0.07, hf =0.4 and R = 92. Stage two lasts until 
about 1.5 MCS, after which S_L{q,t) and Si\{q,t) begin to 
separate and eventually {q, t) decreases towards its stable 
value, (b.) A semi-log plot of S{q,t) for an off-critical quench 
from disorder to clumps with Tf — 0.02, hf — 0.8 and R — 92 
is plotted for comparison. 



Note that in Eq. ([50)) the Qy modes (along the direction 
of the stripes) are decoupled, but the Qx modes (perpen- 
dicular to the direction of the stripes) are not. Hence the 
Fourier modes do not have purely exponential growth as 
in the disordcr-to-order case. 

To solve Eq. (jSOp numerically, we discretize space. For 
a lattice with L x L coarse-grained bins, and therefore 
L X L Fourier modes, we write ^j{qx, Qy) as V-'ij, where Qx 
and Qy are Qx = 2iTi/L and qy = 27rj'/L for Q < i, j < L. 
Because the qy variables are decoupled, we choose a par- 
ticular value j = jo for the remainder of the calculation. 
The discretized Eq. (|50p can be written in matrix form 



dt 



(51) 



where we have left off the j = jo index. The matrix A is 
real and has elements 



A-ik = —MikAk — Sik- 



(52) 



To solve the homogeneous set of equations, dipi{t)/dt — 
Aikipkit), we diagonalize A to find eigenvectors which 
satisfy: 



A a \a a 



(53) 



where the index a — 1,2,...,L. If the eigenvalues A" 
are real and distinct (which we have found to be true for 
the matrices we have numerically diagonalized) then v" 
form a new basis. We expand ip in this basis 



Equation ([5T|) becomes 
d 



dt 



(54) 



(55) 



M{r) — M{rx)- Again, the dissipative dynamics and the 
Glauber dynamics should be equivalent for stage two. 

The nonlinear dynamics of Eq. evolves the initial 
stripe configuration {(jibi'i'x.'f'y^t = 0) = (f>b{rx;t = 0)) 
toward (j)s{fx)- When (j)b{i'xi't) — 4'si'>'x), stage two ki- 
netics begins. In this case stage two kinetics is given by 
Eq. ([35]) where the mobility M{rx) is given by Eq. (gD]) 
with (po{r) — (psi^x) and J = — 1. We find that in order 
to determine the time dependence of the structure factor 
for stage two it is sufficient to consider the dynamical 
equation (|35p without the noise term: 



dt 



~Mirx){A*i/j)ir,t)~^{r,t). 



(49) 



To understand the behavior of the Fourier modes we 
Fourier transform Eq. (|49)) : 



dtpiqx,qy;t) 
dt 



dq'xM{q'x 

i^{qx,qy;t). 



qx)h{q'x,qymq'x,qy]t) 
(50) 



Because u" are linearly independent, Eq. (j55p can be 
written as a set of equations: 



d_ 

dt 



a^it) =a"(OA", 



with solutions 



a"(i) = a"(0)e 



(56) 



(57) 



The structure factor S{qx = 2Tri/L,qy = 27rjo/i;i) = 
^ijoi^) = (l'^ijo(OP) due to the noiseless dynamics is 



Sijo (t) 



V 



E 



a"(0)a" (0) 
V 



where 



— 4'siqx — 2711 /L) and we expect that 
a"(0)a" (0)) ~ V. In the unstable regime, there is 
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at least one positive A" indicating exponential growth. 
From Eq. ([55)) we see the exponential growth lasts 
throughout stage two, but it is growth of eigenvectors 
that will, in general, mix Fourier modes. We expect, 
however, that after some time the term in the sum in 
Eq. ([58]) which corresponds to the largest positive eigen- 
value will dominate, and the growth of Fourier modes will 
be well approximated by pure exponential growth. Note 
that we only expect to see exponential growth for near- 
mean-field systems whose long interaction ranges extend 
the lifetime of the stage two kinetics to r ~ In R at the 
fastest growing eigenvector. 



B. Simulation results 

1. External field quench 

We present the results of Monte Carlo simulations for 
a quench from a stable stripe configuration with hi — 
to ft,/ = 0.8 at Ti = Tf ^ 0.04. Snapshots following this 
quench are shown in Fig. [5] for a system with L — 1024 
and interaction range R ~ 368. We find that the stripes 
quickly become less dense and then narrower. Finally, 
around lOMCS, the stripes begin to form clumps, and 
the modulation along the length of the clumps becomes 
visually perceptible. 

In Fig. [in] we show S{q, t) at various values of q. We 
find that the growth is most rapid along the line of modes 
for which qy = and the line of modes which is closest 
to Qy ~ q* (the peak of the Fourier transform of the in- 
teraction potential.) In Fig. [TO] the Monte Carlo data is 
averaged over at least 700 initial stripe configurations for 
the same parameters as the system shown in Fig. [Hi We 
look at several qx modes with qy q* . There is a steep 
initial increase of S{qx,qy) and then a gradual decrease 
during the first Monte Carlo step, which we interpret as 
stage one kinetics toward the unstable state stripe solu- 
tion. Stage one corresponds to the initial behavior of the 
stripes seen in Fig. [5] followed by narrowing of the stripes. 
Eventually the structure factor modes grow exponentially 
before stage two kinetics breaks down at around 10 MCS. 
Two distinct exponential rates are evident. 

In Fig. [TT] we compare the evolution of the structure 
factor for systems with different interaction ranges. We 
see that the time scale of stage one kinetics is independent 
of the interaction range, and exponential behavior lasts 
longer for larger R, as in the CHC theory. 



2. Understanding stage two kinetics 

To better understand the behavior of the system dur- 
ing stage two, we seek an analogy to the critical quench 
where stage one is entirely absent from the early unstable 
evolution. Thus, instead of quenching the field, we will 
initially prepare our system in the state corresponding to 
the unstable Euler-Lagrange solution. 
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FIG. 9: Several snapshots showing the initial evolution from 
an unstable stripe configuration to a clump configuration in 
the antiferromagnetic Ising model. The lattice size is 1024, 
the interaction range is 368, and Ti = Tf = 0.04. Time is 
measured in Monte Carlo steps beginning after the external 
field quench from hi = to hf = 0.8. The corresponding 
structure factor is shown in Fig. [TOj 
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FIG. 10; Semi-log plot of the structure factor S{q,t) versus 
time in the antiferromagnetic Ising model with R — 368 for 
various values of qx for qy ~ q* following a quench from hi — 
to hf = 0.8 at Ti =Tf = 0.04. The Monte Carlo data is av- 
eraged over about 700 independent stripe configurations. We 
interpret the first Monte Carlo step to be stage one kinetics. 
After this time stage two takes over. At the beginning of stage 
two, the Fourier modes do not grow exponentially. Eventually 
they grow exponentially at two difi'erent rates (starting from 
about 4 Monte Carlo steps). This exponential growth lasts 
until about 10 Monte Carlo steps for R = 368. 
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FIG. 11: Semi-log plot of the structure factor S{q, t) a.t =G 
and qy ~ q* versus time in the antiferromagnetic Ising model 
for interaction ranges 7? = 46, 92, and 368. The quench pa- 
rameters are the same as those in Fig. 1101 We find that the 
time scale of stage one kinetics (about the first MCS) is in- 
dependent of 7?, and the exponential evolution of stage two 
lasts longer for longer R. 



We verify numerically that the unstable stripe solution 
4's{'>'x, fy) to the Euler-Lagrange equation at a particular 
T and h corresponds to the one — dimensional stable 
clump solution (piBirx) with the same temperature and 
external field: (t>s{rx,ry) = (t>s{rx) = 4>i'D{rx)- 

For the stripe to clump transition we used the hf = 
0.8 and Tf — 0.04. We find the stable one-dimensional 
solution by numerically solving Eq. p2p in one dimension 
at this field and temperature. We use the corresponding 
unstable stripe solution 4>s{rx) to find the matrix A in 
Eq. (f52|) which depends on <j>s{rx) through the mobility 

M{r,)=Pf{l-Mrx?)- 

For the parameters considered {hf = 0.8, Tf — 0.04, 
L = 1024, and L/R ~ 2.78), we diagonalize A and find 
two positive eigenvalues of Eq. (Larger lattice sizes 
give a larger number of positive eigenvalues. The depth 
of the quench also effects the number of positive eigen- 
values.) We also find by examining the eigenvectors cor- 
responding to these two positive eigenvalues that each 
Fourier mode has a non-zero contribution from only one 
of the two positive eigenvalues. Thus we expect each 
Fourier mode to be dominated by one of the two expo- 
nential terms with a positive argument and two distinct 
exponential growth rates should eventually be apparent. 
This is consistent with the data shown in Figs.[T^and[T51 

The one-dimensional solution is used to construct ini- 
tial two-dimensional configurations for the Monte Carlo 
simulations. The length of the one-dimensional coarse 
grained lattice is chosen to be the same as the length 
of the two-dimensional Monte Carlo lattice. We use the 
value of 4's{rx) to determine the probability that the cor- 
responding Ising spin is initially up: p — {(j)s{rx) + l)/2. 
We repeat this process over all horizontal strips of the 
Ising lattice. We then run the Monte Carlo simulation 
and monitor S{q, t) for various values of q, averaging over 
independent runs. As expected, the evolution looks simi- 
lar to the evolution of the corresponding S{q, t) in Fig.fTUl 



FIG. 12: Semi- log plot of the structure factor versus time in 
the antiferromagnetic Ising model for various qx values with 
qy — 4.518. The initial configuration is constructed from the 
one- dimensional solution <j!>id to the Euler-Lagrange equation 
at temperature T — 0.04 and field h — 0.8. Data is averaged 
over about 700 independent stripe configurations. Two dis- 
tinct exponential growth rates are apparent from about 3 to 
8 Monte Carlo steps. 
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FIG. 13: A semi-log plot of structure factor for two Fourier 
modes from Fig. [12] along with predictions for linear theory 
exponential growth, fi{t) ~ exp(2Aif) and f2{t) ~ exp(2A2f), 
where the S intercept has been fit to the data. Here Ai and A2 
are the two positive eigenvalues from our matrix calculation. 



with about the first Monte Carlo step removed. In Fig.[T51 
we have plotted two Fourier modes from Fig. [T^ along 
with the exponential curves from Eq. (|58p to indicate 
the growth rates predicted from our matrix calculation. 

After some initial evolution the Fourier modes eventu- 
ally grow exponentially at one of two rates, reflecting the 
coupling of the modes due to the spatial variation of the 
initial state. 

The fact that our theory predicts exponential growth 
rates even though we did not consider the effect of the 
noise on the structure factor indicates that, as in the 
disorder-to-order case, the contribution from the noise 
term has the same time dependence as the noiseless part, 
at least at later times during stage two [3| ■ 
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VII. SUMMARY AND DISCUSSION 

We have presented a generalization of the hnear the- 
ory of Cahn, HiUiard, and Cook [13, [H, IH, [13 for early 
stage unstable evolution. In particular, we have given a 
description of the initial unstable evolution for transitions 
with spatial symmetry breaking and no order parameter 
conservation. One of our predictions is that there are 
two stages of early time kinetics in general for systems 
with long interaction ranges. Stage one is dominated 
by a symmetry preserving nonlinear dynamics and takes 
the system to a symmetry constrained local minimum 
of the free energy on a time scale that is independent 
of the interaction range. Although noise driven fluctua- 
tions may also be growing during stage one, they do not 
evolve significantly in this relatively short-lived regime. 
In stage two these noise driven fluctuations break the 
spatial symmetry of the initial configuration and grow 
for an extended period of time. 

This picture is consistent with the observed behavior 
of the long-range antiferromagnetic Ising model, a sys- 
tem that exhibits several symmetry breaking transitions. 
We find the stage two kinetics of the long-range antifer- 
romagnetic Ising model is characterized by long-lived ex- 
ponential growth of symmetry breaking fluctuations. For 
the disorder-to-order case the symmetry breaking fluc- 
tuations are individual Fourier modes. In contrast, the 
order-to-order case is complicated by the spatial symme- 
try of the initial conflguration, which couples the growth 
of the Fourier modes. For this reason the symmetry 
breaking fluctuations in the order-to-order case are not 
individual Fourier modes, but eigenvectors in Fourier 
space. The evolution of the Fourier modes is described 



by linear combinations of exponential growth and/or de- 
cay, which after some time resembles pure exponential 
growth due to the largest positive eigenvalue. 

We have focused on systems without order parameter 
conservation. A future publication will apply our gen- 
eralized theory for order-to-order transitions with order 
parameter conserving dynamics. 

Our results have several consequences for experimen- 
talists and simulators looking to fit their data to a linear 
theory after an unstable quench. For disorder-to-order 
transitions stage two kinetics, which is analogous to the 
CHC theory, is expected immediately after the quench 
only for systems with conserved order parameters or only 
in special circumstances for systems with no conserved 
order parameters (as in a critical quench.) For order-to- 
order transitions, exponential growth of Fourier modes is 
delayed for two reasons. In stage one exponential growth 
is not expected. In stage two the evolution of Fourier 
modes is represented by a linear combination of expo- 
nential terms, and therefore only approximate exponen- 
tial evolution is expected after some time. Because the 
lifetime of stage two kinetics is limited by the finite inter- 
action range, exponential Fourier mode growth may be 
absent even if the linear evolution of stage two is present. 
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